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Time-of-flight neutron spectroscopy has been used to determine the crystalline electric field (CEF) 
Hamiltonian, eigenvalues and eigenvectors appropriate to the J — 7/2 Yb 3+ ion in the candidate 
quantum spin ice pyrochlore magnet Yb 2 Ti 2 07 . The precise ground state (GS) of this exotic, 
geometrically-frustrated magnet is known to be sensitive to weak disorder associated with the 
growth of single crystals from the melt. Such materials display weak “stuffing” wherein a small 
proportion, 2%, of the non-magnetic Ti 4+ sites are occupied by excess Yb 3+ . We have carried 
out neutron spectroscopic measurements on a stoichiometric powder sample of Yb 2 Ti 207 , as well 
as a crushed single crystal with weak stuffing and an approximate composition of Yb 2 +ccTi 2 -cc 07 + y 
with x — 0.046. All samples display three CEF transitions out of the GS, and the GS doublet itself 
is identified as primarily composed of m j — ±1/2, as expected. However, “stuffing” at low tem¬ 
peratures in Yb 2 +*Ti 2 -* 07 +2/ induces a similar finite CEF lifetime as is induced in stoichiometric 
Yb 2 Ti 207 by elevated temperature. We conclude that an extended strain field exists about each 
local “stuffed” site, which produces a distribution of random CEF environments in the lightly stuffed 
Yb 2 +a;Ti 2 -cc 07 +^, in addition to producing a small fraction of Yb-ions in defective environments 
with grossly different CEF eigenvalues and eigenvectors. 

PACS numbers: 75.25.-j,75.10.Kt,75.40.Gb,71.70.Ch 


I. INTRODUCTION 

Geometrically-frustrated magnetic materials are of 
great current interest due to the diversity of exotic or¬ 
dered and disordered ground states (GSs) that they dis¬ 
play 1 . In particular, cubic pyrochlore magnets with 
chemical composition A 2 B 2 O 7 have been a playground 
for geometric frustration as both the A 3+ and B 4+ sites, 
independently, reside on interpenetrating networks of 
corner-sharing tetrahedra (see Fig. 1), one of the canon¬ 
ical architectures for frustrated ground states in three 
dimension^. The rare earth titanate pyrochlores have 
played a pivotal role in the development of the field, as 
Ti 4+ at the B site is non-magnetic, and the A site can 
be occupied by all trivalent rare earth ions from Sm 3+ to 
Lu 3+ . Many of these rare earth titanates therefore have 
a magnetic A sub-lattice, and the family as a whole gives 
rise to different combinations of magnetic anisotropies 
and interactions, which in turn are responsible for the 
diversity of exotic ground states 3 ! 

As an example, the classical spin ice state has been the 
focus of much attentiorP®, and it results from a combi¬ 
nation of local Ising anisotropy 7 8 and net ferromagnetic 
interaction ^ 8 on the pyrochlore lattice, such that the 
Ising magnetic moments on each tetrahedron obey “ice 
rules” with two spins pointing into each tetrahedron 
and two spins pointing out. analogous to the water ice 
model proposed by Pauling 9 . This results in a sixfold 
degeneracy for a given tetrahedron, and a macroscopic 


degeneracy for the three dimensional network as a whole. 
Yb 2 Ti 207 has attracted much recent attention 10 ^ as 
a candidate for a quantum spin ice ground state | 15 1 16 1 , 
wherein effective S = 1/2 degrees of freedom decorate 
a pyrochlore lattice and interact via net ferromagnetic 
interactions 17 . The GS phase diagram and microscopic 
Hamiltonian appropriate to Yb 2 Ti 2 07 have been ex¬ 
tensively studied^ 19 . The microscopic Hamiltonian 
itself has been determined by modelling spin wave 
dispersion and neutron intensities in the high magnetic 
field, low temperature state of Yb 2 Ti 2 07 ^. This 
work convincingly showed anisotropic exchange to be 
the relevant form of the interactions at low temperatures. 

One of the most interesting features of the GS prop¬ 
erties of Yb 2 Ti 2©7 is its apparent sensitivity to small 
levels of defects that are present in real materials 21 23 . 
Stoichiometric Yb 2 Ti 2 07 is known to display a large and 
sharp anomaly in its heat capacity, C p , near ^ 265 mK. 
However, this anomaly has been observed to be sample 
dependent with samples displaying broader anomalies at 
lower temperatures depending on the exact stoichiometry 
of the material studied EMU! This phenomena is very un¬ 
usual for a three dimensional magnet, as the defect levels 
involved are at the limits of detectability by conventional 
techniques. The variation in stoichiometry, where char¬ 
acterized, is on the order of 1 % level, far removed from 
percolation thresholds in three dimensions. 

Polycrystalline samples of Yb 2 Ti 207 tend to display 
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FIG. 1: The stuffed pyrochlore structure of Yb^Ti^O^. - 

The pyrochlore lattice consists of interpenetrating networks of 
corner-sharing tetrahedra which are generated independently 
by the Yb 3+ (magnetic) and Ti 4+ (non-magnetic) cations. 
The pyrochlore lattice is said to be stuffed wherein a Yb 3+ 
ions, which are normally found on the 16d site (also called 
A-site), also occupy the Ti 4+ 16c site (B-site). 


the sharpest and highest temperature C p anomalies, 
likely due to the lower temperatures required for their 
synthesis, leading to less Ti02 volatization. Single crys¬ 
tals grown by floating zone image furnace techniques typ¬ 
ically display broad low temperature Cp anomalies, or no 
anomalies, and sometimes show multiple peaks in Cp at 
low temperature. Ross et al 21 undertook a detailed neu¬ 
tron crystallographic study of both powder and crushed 
single crystal Yb 2 Ti 2 07 samples, representative of those 
that displayed sharp and broad low temperature Cp sig¬ 
natures, respectively, and showed that the crushed sin¬ 
gle crystal sample displayed weak “stuffing”: a struc¬ 
tural defect in which a slight excess of Yb occupies the 
Ti sub-lattice. Stuffing is illustrated schematically in 
Fig. 1, wherein Yb 3+ ions occupy both the 16d site of 
the Fd3m cubic space group, as well as act as impurities 
on the 16c site normally occupied by Ti 4+ . Weak stuff¬ 
ing was shown to occur at the 2.3% level in the crushed 
single crystals with composition Yb 2 +^Ti 2 -^ 07 +?/ grown 
by floating zone image techniques. In contrast the pow¬ 
der sample grown by solid state synthesis was shown to 
be stoichiometric. 

In this paper we report neutron spectroscopic measure¬ 
ments of the crystalline electric field (CEF) excitations 
in the two powder samples, the stoichiometric powder 
and the crushed single crystal with 2.3 % “stuffing” pre¬ 
viously studied by Ross et al. These measurements allow 
us to accurately determine the eigenvalues and eigenvec¬ 
tors appropriate to the 4 doublets which make up the 
J = 7/2 CEF manifold for Yb 3+ in Yb 2 Ti 2 0 7 . Mea¬ 
surements on both stoichiometric and lightly “stuffed” 
samples allow us to investigate the role of “stuffing” on 
the CEF levels associated with Yb 3+ ions properly resid¬ 


ing on the A-site of the pyrochlore structure. With these 
benchmark measurements in hand, we calculated CEF 
eigenvalues and eigenvectors for the “stuffed” Yb 3+ ions 
residing on the B site, as well as for the A site Yb 3+ ions 
in the presence of oxygen vacancies 


II. CALCULATED CRYSTAL FIELD LEVELS 
FOR Yb 3+ AT THE A-SITE 


Hund’s rules enable the determination of the total an¬ 
gular momentum J of the Yb 3+ ion. The electronic con¬ 
figuration of Yb 3+ is 4/ 13 , resulting in J = | which is 
2 J + 1 = 8-fold degenerate. Within the pyrochlore struc¬ 
ture this degeneracy is lifted by the CEFs at the Yb 3+ 
site due largely to the presence of the eight neighbour¬ 
ing O 2- ions. As illustrated in Fig. 2 (left panel), the 
oxygen environment at the A-site consists of a scaleno- 
hedron, which is a cube distorted along one diagonal that 
forms the local [111] axis. Six oxygen ions, commonly re¬ 
ferred to as 0(2), are located on a plane perpendicular 
to this direction, which is a three-fold rotation axis. The 
other two oxygen ions, referred to as 0(1), are located 
along the local [111] axis in the geometric centre of the 
tetrahedra defined by the A-site Yb 3+ ions. By contrast, 
the environment at the B-site is a trigonal anti prism 
made of six 0(2) oxygen ions surrounding the transition 
metal, as shown in the right panel of Fig. 2. 



A-Site 



B-Site 


FIG. 2: A comparison between the A-site and B-site oxygen 
environment in the pyrochlore structure of Yb^Ti^O^ . - The 
left panel of the figure shows the scalenohedron environment 
generated by the oxygen ions at the A-site where the Yb 3+ 
resides. The symmetry of this structure is similar to that at 
the B-site (right panel) where the Ti 4+ ions are located. 


Following Prather’s convention^, the 3-fold axis 
should be placed along z of the reference system in or¬ 
der to minimize the number of CEF parameters in the 
Hamiltonian. Therefore the resulting CEF Hamiltonian 
for Yb 3+ on the A-site can be written as: 
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Ucef = B%6\ + b\6\ + b\6\ + 

B° 6 d° 6 +Biol + Biol ( 1 ) 

Here we employ the Steven’s operators and CEF 

parameters B™ to approximate the Coulomb potential 
generated by the crystalline electric field due to the neigh¬ 
bouring oxygen atoms. 

The determination of the crystal field parameters B™ 
in Eq. [l] is well suited to inelastic neutron spectroscopy. 
The unpolarized neutron partial differential magnetic 
cross-section can be written within the dipole approxi¬ 
mation as 32 : 


^=c|f(|Q|)S(|Q|,M, (2) 

where U is the scattered solid angle, E' the final neutron 
energy, the ratio of the scattered and incident mo¬ 
mentum of the neutron, C a constant and F(\Q\) is the 
magnetic form factor. The scattering function S(\Q\, Tiuj) 
gives the relative scattered intensity due to transitions 
between different CEF levels. At constant temperature 
and wave vector |Q|, we have: 


S(\Q\,hu) = £ — L(AE + hu>), (3) 

Af e 3 


where a = x, y, z and L(AE ± %uS) = L(Ei — E# ± TmS) 
is a Lorentzian function which ensures energy conserva¬ 
tion as the neutron induces transitions between the CEF 
levels i i', which possess a finite energy width or life¬ 
time. 

The procedure for fitting our inelastic neutron scatter¬ 
ing (INS) data assumes an initial set of CEF Hamiltonian 
parameters. The Hamiltonian is then diagonalized to find 
the corresponding CEF eigenfunctions and eigenvalues. 
S(\Q\, Tiuj) is then computed and directly compared with 
the experimental results and, finally, the CEF parame¬ 
ters are tuned such that the y 2 between the calculated 
and measured S(\Q\,hw) is minimized. 


III. EXPERIMENTAL DETAILS 

Now we turn our attention to the details of the INS 
experiment which was conducted on two samples of 
Yb 2 Ti 207 of different stoichiometry. We will refer to 
these two samples as the stoichiometric (x = 0) and 
stuffed powder (x = 0.046) samples hereafter. Specific 
details regarding their synthesis and characterization can 
be found in previous work by Ross et al. 21 . These sam¬ 
ples were studied utilizing the SEQUOIA direct geome¬ 
try time-of-flight spectromete^, which is located at the 
Spallation Neutron Source at Oak Ridge National Lab¬ 
oratory. In an INS experiment the CEF excitations are 


manifested as dispersion-less features with the strongest 
scattering intensity expected at low \Q\ positions as a 
result of their magnetic origin. SEQUOIA is the ideal 
instrument for the investigation of the CEFs since it 
provides low \Q\ coverage and a large, dynamic (|Q|,E) 
range. 

12g of each of the stoichiometric and stuffed powder 
samples were loaded in an aluminium flat plate with di¬ 
mensions 50mm x 50mm x 1mm and sealed with indium 
in helium atmosphere. An empty can with the same di¬ 
mensions was loaded with the two samples on a three 
sample changer in a closed-cycle refrigerator. Measure¬ 
ments have been performed over a range of tempera¬ 
tures from T = 5 K to 300 K, and utilizing neutrons 
with incident energy Ei = 150meV giving an elastic en¬ 
ergy resolution of ±2.8 meV. The energy resolution im¬ 
proves with increasing energy transfer and is ~1.4meV 
and ^lmeV for energy transfers of 80 and 115meV, re¬ 
spectively. The corresponding chopper settings selected 
were To = 120Hz and FCi = 600Hz. Similar measure¬ 
ments were conducted on the empty can for use as a 
background measurement. 


IV. INELASTIC NEUTRON SCATTERING 
FROM CRYSTAL FIELD EXCITATIONS IN 

Yb2Ti207 


The INS spectra taken on both the stoichiometric and 
stuffed powder samples are shown in Fig. 3 for energies 
up to Ei = 150meV at T = 5 K. An empty can data set, 
taken as background at the same temperature, has been 
subtracted from both the sample data sets. As seen in 
Fig. 3, we observe a set of dispersion-less excitations 
arising from both the CEF transitions, as well as optical 
phonons. Below 60 meV, the intensity associated with 
the dispersion-less excitations increases with increasing 
| Q |, characteristic of inelastic scattering from phonons. 
In contrast, the excitations at 76.7 meV, 81.8 meV and 
116.2 meV, highlighted with blue, grey and red arrows, 
respectively, in Fig. 3, show inelastic scattering which 
increases with decreasing |Q|, consistent with magnetic 
scattering. We therefore ascribe these three dispersion¬ 
less excitations with CEF transitions from the GS dou¬ 
blet. 

Yb 3 + possesses thirteen electrons in its almost filled 
4/ shell and, as a consequence of Kramers’ theorem, its 
eight-fold degenerate CEF levels can be maximally split 
into four doublets. We associate the three magnetic ex¬ 
citations observed at T = 5 K in both the stoichiometric 
and the stuffed powder samples shown in Fig. 3, with the 
transitions between the CEF GS doublet and the CEF 
excited state doublets, as indicated in the inset of Fig. 4. 
These three transitions between CEF doublets account 
for all the CEF states within this J = 7/2 multiplet ap¬ 
propriate to Yb 3+ . As will be described in further detail, 
the CEF transitions observed at low temperatures in the 
stoichiometric sample are sharper in energy than those 
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FIG. 3: Inelastic neutron scattering spectra- S(\Q\,hu) ob¬ 
tained for the stoichiometric powder and the stuffed powder 
samples at T = 5 K are shown in panels (a) and (b) respec¬ 
tively, with the corresponding T — 5 K empty can subtracted 
from each data set. The three horizontal arrows in blue, gray 
and red highlight the three crystal field excitations which are 
found at 76.7 meV, 81.8 meV and 116.2 meV, respectively. 


in the presence of light stuffing. We shall restrict our 
quantitative analysis of the CEF spectra to the case of 
the stoichiometric powder sample. Therefore we use the 
x = 0 data set in Fig. 3 a), and perform a \Q\ integrated 
cut (|Q| = [4.5,5.25]A _1 ), yielding the neutron scatter¬ 
ing intensity as a function of energy. This data set is 
shown in the main panel of Fig. 4, wherein the intensity 
at the peak of the CEF transition at 81.8 meV has been 
normalized to unity. The relative intensities of the CEF 
transitions at 76.7 meV, 81.8 meV and 116.2 meV, as well 


as the energy of these transitions from the GS, constrain 
the CEF Hamiltonian. 

The cut shown in Fig. 4 has been fitted to a model 
for the inelastic scattering 32 ' arising from dipole allowed 
transitions between the CEF GS doublet and the three 
excited states, using Eq. 3. The starting parameters in 
the CEF Hamiltonian for Yb 3+ were those determined by 
Bertin et a/P^for Yb 2 Ti 2 C >7 within the point charge ap¬ 
proximation. This calculation yielded CEF transitions at 
^60 meV, ^70 meV and ~90 meV, in contrast to those 
determined experimentally in Fig. 3. These starting pa¬ 
rameters were then refined and a best fit (shown in solid 
red) to the \Q\ integrated cut in Fig. 4 was obtained. 
All three CEF transitions were fit using the same reso¬ 
lution determined energy width. In addition, relatively 
weak Lorenzian lineshapes near 56 meV, 70 meV and 
100 meV phenomenologically describe the inelastic scat¬ 
tering from the phonons which are in near proximity to 
the CEFs, and improved the fit when included. As can 
be seen, the overall description of the INS data from the 
stoichiometric sample in Fig. 4 is very good. 



0.0 — 1 — a — 1 —■ 1 ■— 1 —*— 1 —*— 1 — 

60 70 80 90 100 110 120 
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FIG. 4: |Q| integrated cut for the stoichiometric powder at T 
= 5 K. -A \Q\ integrated cut (|Q| = [4.5,5.25]A -1 ) obtained 
from the INS spectra in Fig.3(a) for the stoichiometric powder 
at T — 5 K is shown. The arrows indicate the corresponding 
positions of the CEFs. The inset outlines the corresponding 
CEF transitions from the GS doublet. 

The values which were obtained for the CEF Hamil¬ 
tonian parameters from the fit relevant to Yb 3+ at the 
A-site are given in Table [I] and the resulting energy eigen¬ 
values and eigenvectors are given in Table [H| The GS 
doublet for Yb 3+ is comprised primarily of m j = db 1/2 
and the corresponding low temperature anisotropic g- 
tensor components are given by g±_ = 3.62 db 0.15 and 
g z = 1.85 ± 0.10, where z corresponds to the local [111] 
axis. The error bars on the g values were estimated by 
exploring the sensitivity of the fitting procedure on the 
low temperature stoichiometric powder data, to differ- 
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B™ (meV) 

Calculated 

Fitted 

Ratio 

B u 2 

1.270 

1.135 

0.894 

B° a 

-0.0372 

-0.0615 

1.653 

B\ 

0.275 

0.315 

1.145 

Bl 

0.00025 

0.0011 

4.4 

Bl 

0.0023 

0.037 

16.087 

B% 

0.0024 

0.005 

2.083 


TABLE I: A comparison of the calculated crystal field pa¬ 
rameters ( B with those obtained by fitting INS data from 
the stoichiometric powder at T = 5 K. 


ent quantitative descriptions of the background phonons. 
Note that this determination of the g-tensor is performed 
at zero magnetic field, in contrast to, for example, the 
analysis of spin wave data in high magnetic fields 20 ! The 
g-tensor is consistent with previous estimates for Yb 3+ 
in Yb 2 Ti 2 07 ®^, and the dominant mj = db 1/2 char¬ 
acter of the GS validates the effective S = 1/2 quantum 
description for the Yb 3+ moment. In order of ascending 
energy, the excited state doublets correspond primarily 
to mj = db 7/2, mj = db 3/2, and rrij = db 5/2. 



FIG. 5: Energy cuts of the \Q\ = [4.5,5.25]A~ x integrated 
inelastic scattering for the stoichiometric powder sample as a 
function of temperature. - Energy cuts were taken with an 
appropriate background subtraction and all energy cuts for 
temperature above T — 5 K has been vertically translated for 
clarity. 


A. Temperature and weak stuffing dependence of 
the A-site CEF transitions 


Inelastic neutron scattering measurements have also 
been carried out on both the stoichiometric and stuffed 
powder samples as a function of temperature, using 
Ei = 150meV neutrons. Energy scans of the \Q\ = 
[4.5,5.25]A _1 integrated inelastic scattering for the sto¬ 
ichiometric powder sample are shown as a function of 
temperature in Fig. 5. The T = 5 K data set is the same 
x = 0 data set shown in Fig. 4, and the peak intensity 
associated with the 81.8 meV CEF transition at T = 5 
K has been normalized to unity. 

We note that the maximum temperature employed in 
these measurements, 300 K, corresponds ~ 27 meV and 
consequently only the GS doublet of Yb 3+ is substan¬ 
tially occupied at any temperature. Qualitatively, the 
thermal fluctuations have three effects on the CEF neu¬ 
tron spectra: the CEF excitations broaden appreciably 
in energy; the maximum peak intensity diminishes; the 
energy of the CEF excitations softens slightly. This lat¬ 
ter effect is somewhat subtle, but it can be seen in Fig. 5 
by drawing a fiducial dashed line positioned at the centre 
of the CEF transitions at T = 5 K. 

The energy width of the CEF excitations can be quan¬ 
titatively examined by fitting the data sets as in Fig. 4, 
but now with a damped harmonic oscillator (DHO) line- 
shape for the three CEF transitions, each with the same 
temperature dependent energy width. At the energy 
transfers and temperatures of interest, the DHO can be 
approximated by a single Lorentzian for each mode.The 
form of this line-shape is given by: 


L(E) = 1 (^r) 

V ; n(E- A£) 2 + (^*) 2 ’ 

which is a Lorentzian function of energy with width r o bs 
and centered on the energy of the CEF transition A E. 
The falloff of the maximum peak intensity as a function 
of increasing temperature along with the concomitant 
broadening in the energy widths of the transitions en¬ 
sures that the integrated spectral weight of this inelastic 
scattering is almost temperature independent, consistent 
with the CEF transitions being from the GS doublet of 
Yb 3+ at all temperatures considered here. The common 
energy width or inverse lifetime of the three CEF exci¬ 
tations extracted from this analysis is plotted as a func¬ 
tion of temperature in Fig. 6, where we removed the 
resolution contribution to the widths using the following 
relation : 


rfntrinsic(T) = ^(T) - (5) 

Figure 6 then shows pronounced growth in the intrin¬ 
sic energy width of the CEF transitions from ^3.5 meV 
at T=5 K to ^ 9 meV at T=300 K. This growth mirrors 
the temperature dependence of the Yb 3+ mean squared 
displacements (MSDs) in stoichiometric Yb 2 Ti 2 07 as de¬ 
termined by powder neutron diffractiori 21 , and this is also 
reproduced in Fig. 6 for direct comparison. It is worth 
noting that we plot one of the anisotropic components 
of the atomic displacement parameter, Un, which repre¬ 
sents the anisotropic MSD of the Yb 3+ ions from their 
average positions. Although, for the specific case of the 
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TABLE II: The CEF eigenvalues and eigenvectors for Yb 3+ at the A-site of Yl^TLOr. The first column displays the CEF 
eigenvalues of the system, while the corresponding eigenvectors are given in each row in terms of the mj basis. 


y b 3 + 

ion in the stoichiometric powder, it can be shown 
that Un is equal to the other components of the atomic 
displacement, so that the MSD is effectively isotropic. 
The inset to Fig. 6 shows the slight softening of the CEF 
energies with increasing temperature. The softening is 
most pronounced for T < 125 K, and then it flattens out 
at higher temperature. The significance of 125 K = 11.3 
meV is not completely clear, although it could be related 
to the temperature scale associated with the top of the 
acoustic phonon band in Yb 2 Ti 207 , which is typically 10 
-15 meV for transition metal oxides. This effect is rela¬ 
tively small with a maximum observed softening of ^ 0.8 
meV at T = 300 K compared with T=5 K. 



Temperature (K) 

FIG. 6: Energy width and shift of the CEF transitions as 
a function of temperature- The intrinsic energy width of the 
CEF transitions as extracted using Eq. 5. The intrinsic en¬ 
ergy width corresponding to the stuffed powder is higher at 
both 5 K and 300 K in comparison to the stoichiometric sam¬ 
ple. The green dots shows the mean square displacement 
(MSD) of the Yb 3+ taken from rel 21 ^and scaled in such a way 
that the MSDs and the energy widths of the CEF excitations 
in the stoichiometric powder have the same value at T=100K. 
The inset shows the shift in energy of the CEF transitions as 
a function of temperature, relative to that at T— 5 K. 

The finite intrinsic energy widths of the CEF excita¬ 
tions in stoichiometric powder must originate from dy¬ 


namics in the lattice not captured by the static structure. 
The strong resemblance of the temperature dependence 
of the Yb 3+ MSDs to that of the CEF energy widths 
suggests that zero point fluctuations determine the low 
temperature energy widths, and the thermal population 
of phonons give rise to the larger widths at finite temper¬ 
ature. The ionic displacements associated with both zero 
point fluctuations and the phonons will distort the local 
environment at the Yb 3+ site, giving rise to a distribution 
of CEF transitions energies. Displacements of harmonic 
phonons time-average to zero, consequently the distribu¬ 
tion of CEF transitions in the presence of the phonons 
is approximately centred on the zero temperature CEF 
transitions, and the transitions remain well defined at all 
temperatures. 

Figure 7 shows a direct comparison of the T = 5 K, 
|Q | = [4.50,5.25] A -1 integrated inelastic scattering for 
the stoichiometric and stuffed samples. This comparison 
clearly shows the CEF excitations in the lightly-stuffed 
sample have a considerably larger energy width even at 
low temperatures. The INS from the stuffed powder was 
also fit to the same DHO line-shape as was fit to the sto¬ 
ichiometric data, and the corresponding intrinsic energy 
widths at T = 5 K and 300 K are also plotted in Fig. 
6. The intrinsic energy widths of the CEF excitations in 
the stoichiometric and lightly stuffed samples differ by 
only ^ 5 % at T = 300 K; this is not surprising as the 
intrinsic energy widths in either sample are expected to 
be dominated by thermal fluctuations. However, at T = 
5 K, the intrinsic energy width of the stuffed sample is ^ 
5.2 meV - approximately that displayed by the stoichio¬ 
metric sample at ~ 125 K. 

The finite energy widths of the CEF excitations at low 
temperature in the lightly stuffed sample demonstrate 
that weak stuffing must induce a relatively large “vol¬ 
ume of influence” about each “stuffed” Yb ion, that is, 
about each Yb ion residing on a B-site. The correspond¬ 
ing distortion field about each stuffed site gives rise to a 
distribution of A-site Yb 3+ environment, similar to those 
associated with phonons at finite temperature. Indeed, 
the stuffed B-site Yb ions would also give rise to CEF 
spectra completely distinct from those of Yb 3+ in A-site 
environments. However the concentration of such stuffed 
B-site Yb ions is known to be small, ~ 2.3 %, in our 


























7 



FIG. 7: A direct comparison between the INS from CEF 
transitions in the stuffed powder (in orange) and the stoichio¬ 
metric powder (in green) Yb 2 Ti 2 07 samples. - Energy cuts 
taken at T = 5 K within a range \Q\ = [4.50,5.25]A _1 for the 
stuffed powder (orange) and stoichiometric powder (green) 
samples. 


lightly stuffed sample and we were unable to find con¬ 
vincing evidence for such excitations in this study. 


T = 5 K is shown in Table [TTI Tables III and IV show 
the CEF eigenvectors and eigenvalues for Yb ions in two 
different impurity sites: one “stuffed” on the B-site and 
one on the A-site but in the presence of one 0(1) oxy¬ 
gen vacancy. Both of these impurities are expected to be 
present at some small concentration in the lightly stuffed 
sample. In principle 0(2) vacancies can also occur, but 
these are higher energy defects than 0(1) vacancies^. 
These CEF eigenvectors and eigenvalues have been cal¬ 
culated within the point charge approximation and then 
corrected with the same ratio of B™ parameters as were 
determined from the comparison between the INS fit¬ 
ted values and the calculated values in stoichiometric 
Yb 2 Ti 207 at low temperatures. 

The differences between the fitted CEF parameters and 
the calculated ones originate mainly from the overlap of 
the 4f orbital of the rare earth with the 2p orbitals of 
the ligands. The corresponding effects on the B™ terms 
are difficult to estimate from first principles and are not 
considered here. Instead we calculate the CEF param¬ 
eters for the defective environments, and then correct 
these with the same B™ ratios which were determined 
in the treatment of CEFs for the stoichiometric sample. 
Finally, we verified that this correction did not dramati¬ 
cally affect the properties of the system; that is that the 
CEF properties within the defective environments did 
not qualitatively depend on this correction. 


V. CRYSTAL FIELD CALCULATIONS FROM 
FIRST PRINCIPLES WITHIN THE POINT 
CHARGE APPROXIMATION 

As discussed in Sec. [IlJ the CEF Hamiltonian param¬ 
eters can be fitted using INS data. In principle, these 
parameters can also be calculated from first principles 
within a point charge model following the procedure out¬ 
lined in Refs.OIMSSrhe CEF interactions are treated as a 
perturbation to the spin-orbit coupling and the resulting 
eigenfunctions of the CEF Hamiltonian are expressed as 
a linear combination of the | J, mj) states within the 2 i± 
manifold multiplets. 

We have calculated the CEF Hamiltonian parameters 
in this manner and the ratio between these ab initio pa¬ 
rameters and those extracted from fit of the INS data 
are shown in Table [T| Note that the ionic positions were 
taken from the crystallographic refinement of stoichio¬ 
metric Yb 2 Ti 2 07 21 . As can be seen, apart from the Bq 
and Bq terms which are relatively small, all the CEF pa¬ 
rameters agree well with those extracted from fitting the 
INS data. This result is remarkable given the simplicity 
of the point charge approximation, and it gives us con¬ 
fidence that we can make reasonable predictions on the 
strength of the CEF Hamiltonian parameters determined 
in this way. 

The CEF eigenvectors and eigenvalues determined 
from the fit to the INS in the stoichiometric powder at 


A. Crystal Field Calculation for Yb 3+ at the A-site 
in an Oxygen Depleted Environment 


We first consider an A-site Yb 3+ in the presence of a 
single 0(1) oxygen vacancy. This depleted environment 
is expected to break the symmetries of the CEF Hamil¬ 
tonian. However if Prather’s convention is satisfied, only 
the inversion operation is lost and the number of CEF 
parameters in our Hamiltonian is unchanged from the 
stoichiometric case. Note that in general, a broken sym¬ 
metry can dramatically affect the CEF Hamiltonian, and 
it may be necessary to add terms to the CEF Hamiltonian 
to better approximate the defective Coulomb potential. 

The eigenvalues and eigenvectors calculated for A-site 
Yb 3+ in the presence of a single 0(1) oxygen vacancy 
are shown in Table [TTT| As previously discussed, we have 
performed the point charge calculation and then scaled 
these results by the B™ ratios taken from the fitted and 
calculated CEF Hamiltonian terms determined for the 
stoichiometric powder. 

Kramers’ degeneracy still protects the Yb 3+ ion in this 
defective environment and all the CEF eigenvectors ap¬ 
pear as doublets. However, the gap between the GS and 
the first excited state is much larger than in the stoichio¬ 
metric case. Moreover an examination of the eigenvectors 


in Table III shows that the GS is now a pure linear com¬ 
bination of mj = ±7/2 states. These GS eigenvectors are 
similar to what is found in stoichiometric Dy 2 Ti 207 crys¬ 
tals, where the GS is comprised of a linear combination 
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TABLE III: The CEF eigenvalues and eigenvectors calculated for Yb 3+ in a depleted oxygen environment at the A-site in 
Yb 2 Ti 207 - The first column displays the CEF eigenvalues of the system, while the corresponding eigenvectors are given in 
each row in terms of the mj basis. The ratio correction was applied in order to arrive at the final CEF parameters. 


of the maximal mj states, pure raj=±15/2. 



FIG. 8: Magnetization of the A-site Yb 3+ ion in presence 
of an 0(1) oxygen vacancy. - A coloured sphere represents 
the tip of the magnetization vector centred at the Yb 3+ ion 
in response to an applied external field of H = 1 T which 
rotates within the plane normal to [111]. Many values of 
the magnetization are shown simultaneously for applied fields 
uniformly distributed on the unit sphere. For the same [111] 
component of magnetization, the centre of the spheres are 
sufficiently close to each other, that the final impression is 
that of rings. The tightness of the rings around the [111] axis 
implies stronger Ising-like behaviour of the Yb 3+ moment, in 
this case with a calculated magnetic moment of 3.953 jiB along 
the local [111] direction. 


With the CEF eigenfunctions in hand, we can calcu¬ 
late the local magnetization of the A-site Yb 3+ in the 
presence of a single 0(1) oxygen vacancy. This is shown 
in Fig. 8 where the magnetization anisotropy for Yb 3+ 


moment within the defective scalenohedron is displayed. 
In both Figs. 8 and 9, coloured spheres represent the tips 
of magnetization vector centred at the Yb 3+ ion, in re¬ 
sponse to a given applied external field of H = 1 T with a 
particular [111]-component and which precesses around 
the plane normal to [111]. Many values of the magne¬ 
tization are shown simultaneously for applied fields uni¬ 
formly distributed on the unit sphere. The colour scale 
is chosen in order to match the strength of the magnetic 
moment from 0/i# (dark blue) to 4 hb (dark red). For the 
same [111] component of magnetization, the centre of the 
spheres are sufficiently close to each other, that the final 
impression is that of rings (this is seen more clearly in Fig. 
9). The tightness of the rings about the [111] axis implies 
stronger Ising-like behaviour of the Yb 3+ moment. For 
the case of the A-site Yb 3+ in the presence of a single 
0(1) oxygen vacancy in Fig. 8, the calculated magnetic 
moment along the local [111] direction is 3.953/ig, close 
to the full moment of 4/i#. This suggests negligible pre¬ 
cession around this easy axis and thus strong Ising-like 
behavior. The local magnetization anisotropy for the A- 
site Yb 3+ in the presence of a single 0(1) oxygen va¬ 
cancy is dramatically different fro m the st oichiometric 
case, where it was planar or XY-likd ^ 6 * 37 * 41 ( and is sim¬ 
ilar to that displayed by the classical spin ice magnets 
Dy 2 Ti 2 07 and Ho 2 Ti 2 07 . This anisotropy is present even 
without the B™ ratio corrections and thus it depends 
only on the depleted environment. 


B. Crystal Field Calculation for “Stuffed” Yb 3+ at 
the B-site 


Neutron powder diffraction measurements on the 
stuffed powder samples, derived from crushed single 
crystals, show such Yb 2 Ti 2 07 samples to be lightly 
stuffed at the 2.3 % level. We therefore calculate the 
CEF eigenvectors and eigenvalues for “stuffed” Yb 3+ 
at the B-site. This is an interesting case to consider 
since the local environment at the B-site is very different 
from that considered to this point for A-site Yb3±. 
Nonetheless, the symmetry at the B-site is similar to the 
A-site, and we can simply rotate the trigonal anti-prism 
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cage around the B-site so that the local [111] direction is 
aligned along z and the C 2 axis along y. Therefore the 
form of the B-site Yb 3+ CEF Hamiltonian is identical to 
Eq. [l] The calculated CEF eigenvectors and eigenvalues 
for B-site Yb 3+ are displayed in Table IV Once again 
we have performed the point charge calculation and then 
scaled these results by the B™ ratios determined for the 
stoichiometric powder. 


Kramers’ degeneracy is again invoked and all CEF lev¬ 
els are again doublets, now with a gap of 105 meV be¬ 
tween the GS and the first excited state. The eigen¬ 
vectors within the GS doublet are again almost pure 
mj=±7/2 spin states, and we expect the anisotropy to 
be Ising like. This is borne out by again applying a small 
rotating external magnetic field to probe the shape of 
the anisotropy of the GS magnetic moment as was done 
for Yb 3+ at the A-site in an oxygen depleted environ¬ 
ment. Following the previous convention, Fig. 9 shows 
spheres centred along the direction in which the magnetic 
moment points. The colour scale matches the magnetic 
moment size in units of /i#. In contrast with the previ¬ 
ous case, the spheres are not located exactly along the 
local [111] direction, thus we anticipate a small precession 
about the easy axis with the maximum value of the mag¬ 
netic moment given by fi = 3.25/r#- This is similar to 
what was calculated for the defective scalenohedron as¬ 
sociated with Yb 3+ at the A-site in an oxygen depleted 
environment, but the anisotropy is not as Ising-like as it 
was for that case. 

It is interesting that the influence of oxygen vacancies 
on A-site Yb 3+ ions and that of stuffed B-site Yb 3+ ions 
have similar effects on the CEF eigenvalues and eigen¬ 
vectors. Stuffing Yb 3+ onto B-sites nominally occupied 
by Ti 4+ will require oxygen vacancies to preserve charge 
neutrality, and we expect both phenomena to occur at 
some level in single crystals grown from the melt. The 
symmetry of the defective scalenohedron, especially when 
two 0(1) are missing, is reasonably similar and to the 
trigonal anti-prism associated with stuffed Yb 3+ at the 
B-site, so it is perhaps not so surprising that these envi¬ 
ronments produce similar effects on the CEF eigenvalues 
and eigenvectors of the magnetic ions. On the other hand 
the similarity of the CEFs associated with the defective 
sites may complicate the analysis of experimental data 
from sufficiently defective non-stoichiometric Yb 2 Ti 207 , 
as the measured spectrum may show overlap originating 
from at least two distinct defective environments. 


VI. CONCLUSIONS 

Neutron scattering measurements on the stoichiomet¬ 
ric and stuffed powder samples of the quantum spin ice 
candidate system Yb 2 Ti 2 0 y have been carried out to 
probe the CEF eigenvalues and eigenvectors associated 
with the J = 7/2 Yb 3+ ion in these environments. Anal¬ 
ysis of the INS from the stoichiometric powder show the 



FIG. 9: Magnetization of the stuffed” Yh 3+ ion at the 13- 
site of Yb 2 Ti 207 . - Each sphere represents the tip of the 
magnetization vector centred at the Yb 3 + ion in response to 
an applied rotating external field of H = 1 T. Neighboring 
contiguous spheres form ring pattern. The ellipsoidal shape 
of the ring pattern suggests Ising-like behaviour of the rare 
earth moment along the local [111] direction, however the 
transverse extent of the ellipsoidal shape indicates a small 
precession around the Ising axis. The calculated magnetic 
moment is 3.25 /is- 


GS doublet to correspond primarily to mj = ± 1/2, with 
anisotropic g-tensor components given by g±_ = 3.62 ± 
0.15 and g z = 1.85 db 0.10, which therefore imply local 
XY anisotropy as expected. The eigenvalues and eigen¬ 
vectors of the three excited state CEF doublets were all 
identified above a lowest CEF gap of ^ 76.7 meV at T 
= 5 K. The energy widths, or inverse lifetimes of the 
CEF states broaden with increasing temperature from ^ 
3.5 meV at low temperatures up to the highest measured 
temperature, 300 K. This broadening mirrors the temper- 
arure dependence of the Yb 3+ MSDs, and we associate 
it with the effects of zero point motion and phonons on 
the CEF environment around the A-site Yb 3+ site. In¬ 
terestingly, the corresponding measurements and analy¬ 
sis of the lightly stuffed powder sample (x = 0.046) of 
Yb 2 Ti 2 07 shows a similar set of CEF transitions, how¬ 
ever the energy width or inverse lifetime of the transi¬ 
tions are intrinsically further broadened even at T = 5 
K. This effect is ascribed to a relatively large strain field, 
or “volume-of-influence”, associated with each stuffed B- 
site Yb 3+ . 

We have also calculated the CEF eigenvectors and 
eigenvalues associated with Yb 3+ in two defective envi- 
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TABLE IV: The CEF eigenvalues and eigenvectors calculated for Yb 2 Ti 207 at B-site. The first column displays the Crystal 
Field spectrum of the system, while the corresponding eigenvectors are given in each row in terms of mj basis. The ratio 
correction was applied in this case to the CEF parameters. 


ronments of relevance for non-stoichiometric Yb 2 Ti 207 . 
These are for Yb 3+ at the A-site in an oxygen depleted 
environment, and for stuffed Yb 3+ at the B-site. These 
calculations were performed within the point charge ap¬ 
proximation and make use of a comparison of the cal¬ 
culated and measured eigenvalues and eigenfunctions in 
the stoichiometric sample to benchmark these results. 
Both of these defective environments give rise to related 
effects wherein the GS doublet is now primarily made 
up of mj = ± 7/2, and the local anisotropy of the GS 
eigenfunctions are Ising-like, although with stronger Ising 
anisotropy for the oxygen depleted environment than for 
case of the stuffed Yb 3+ at the B-site. 

These results put the nature of the CEF eigenfunc¬ 
tions and eigenvalues associated with Yb 3+ in the quan¬ 
tum spin ice candidate system on a much firmer foot¬ 
ing and indicate how the nature of the GS moment 
and anisotropy in this system can be sensitive to small 


amounts of disorder, as is known to occur in single crys¬ 
tal samples grown from the melt using floating zone tech¬ 
niques. 
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